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ABSTRACT 

We present VLA H i observations and UH88" deep optical B- and i?-band observations of the IR 
luminous merger Arp 299 (= NGC 3690 + IC 694). These data reveal a gas-rich (Mhi = 3.3 x IO^Mq) 
optically faint (/x_b^27 mag arcsec"^, ^^^26 mag arcsec"^) tidal tail with a length of over 180 kpc. 
The size of this tidal feature necessitates an old interaction age for the merger (^750 Myr since first 
periapse), which is currently experiencing a very young star burst (;^20 Myr). The observations reveal 
a most remarkable structure within the tidal tail: it appears to be composed of two parallel filaments 
separated by ~20 kpc. One of the filaments is gas rich with little if any starlight, while the other is 
gas poor. We believe that this bifurcation results from a warped disk in one of the progenitors. The 
quantities and kinematics of the tidal H i suggest that Arp 299 results from the collision of a retrograde 
Sab-Sb galaxy (IC 694) and a prograde Sbc-Sc galaxy (NGC 3690) that occurred 750 Myr ago and which 
will merge into a single object in ^ 60 Myr. We suggest that the present IR luminous phase in this 
system is due in part to the retrograde spin of IC 694. Finally, we discuss the apparent lack of tidal 
dwarf galaxies within the tail. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 3690, IC 694, Arp 299) — galaxies: 
interactions — galaxies: ISM — galaxies: kinematics and dynamics — galaxies: 
peculiar — galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 

Arp 299 is a nearby {Vhei—SO&O km s~-^)^ pecuhar sys- 
tem comprised of two highly distorted disk galaxies, with 
NGC 3690 (= Mrk 171A = UGC 6471 = VV118b) to the 
west and IC 694 (= Mrk 171B = UGC 6472 = VV118a)3 to 
the east. There is also a compact spheroidal galaxy lying 
approximately 1' to the NW which is a the same velocity 
as Arp 299 (MCG-hlO-17-2a^ Vhe/=3100 km s'^, FairaU 
1971, Sargent 1972). Figure [| presents a i?-band image of 
this system, taken from the optical data described later in 
this paper. Each of the above named systems are labeled, 
as is the unrelated background bridge-tail-bar system Arp 
296. 

Arp 299 appears to be in an advanced stage of merging, 
with the two disks in contact, but the nuclei still separated 
by about 20" (4.7 kpc)''. This impression is supported by 
radio continuum and near-infrared (NIR) imaging (Gehrz, 
Sramek & Weedman 1983) and millimeter spectral line 
observations (Casoli et al. 1989, 1999; Sargent & Scoville 
1991, Aalto et al. 1997), which locate two major concen- 
trations of molecular gas at peaks in the radio and NIR 
(regions A and B in Fig. |^, using the naming convention 
of Gchrz, Sramek & Weedman 1983). Source A is identi- 
fied with the nucleus of IC 694, while source B is identified 



with the nucleus of NGC 3690.^ Another concentration of 
molecular gas occurs at the region of disk overlap (labeled 
C-C in Fig. |l|). 

Spectroscopic observations of the NIR peaks in radio re- 
combination lines (Anantharamaiah et al. 1993; Zhao et al. 
1997), the Mid infrared (Dudley & Wynn-Wilhams 1993; 
Dudley 1998), the NIR (Gehrz, Sramek & Weedman 1983; 
Beck, Turner & Ho 1986; Nakagawa et al. 1989; Ridge- 
way, Wynn-Wilhams & Beckhn 1994; Doherty et al. 1995; 
Smith et al. 1996; Langon, Rocca-Volmerang & Thuan 

1996) , optical recombination lines (Bushouse & Gallagher 
1984; Keel et al. 1985; Armus, Heckman & Miley 1989) 
and in the UV (Augarde & Lequeux 1985; Kinney et al. 
1993; Roberts 1996; Vacca et al. in preparation) fail to 
reveal any definitive evidence for an AGN, and are con- 
sistent with the dominant energy source being a compact 
starburst (although there is evidence for some contribu- 
tion due to a compact source; see Shier, Rieke & Rieke 
1996; Carral, Turner & Ho 1990; Jones et al. 1990; Jones 

1997) . This picture is consistent with continuum observa- 
tions in the radio (Condon et al. 1982; Gehrz, Sramek & 
Weedman 1983; Condon et al. 1990, 1991; Smith, Lons- 
dale & Lonsdale 1998), far infrared (Joy et al. 1989), NIR 
(Telesco, Decher & Gatley 1985; Carico et al. 1990; Wynn- 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

^AU radial velocities quoted in this paper are heliocentric. 

^There is some confusion as to the naming of this system, wrhich is addressed separately in the appendix. In the following, we following the 
naming convention used most widely in the literature. 

^Adopting the distance of 48 Mpc from Sanders, Scoville & Soifer (1991), which assumes Ho=75 km Mpc~^ and the Virgocentric flow 
model of Aaronson et al. (1982). At this distance l'=14 kpc. 

^High resolution NIR imaging by Wynn- Williams et al. (1991) shows that source B is itself a double, separated by about 3". 
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Willams et al. 1991; Miles et al. 1996), and X-ray (Rieke 
1988, Zezas et al. 1998). 

Within the main body, vigorous star formation is taking 
place, with an inferred star formation rate (SFR) of ^^50 
M0yr~^ (Heckman, Armus & Miley 1990), and an age of 
^20 Myr (Augarde & Lequeux 1985; Stanford 1989; Naka- 
gawa et al. 1989; Meurer et al. 1995; Vacca et al. in prepa- 
ration). A deep Ha -I- [Nil] image of this system reveals 
a "magnificent large-scale filamentary structure made up 
of arcs and loops that surround the galaxies at radii ~ 6- 
12 kpc" (Armus, Heckman & Miley 1990). Recent x-ray 
observations reported by Heckman et al. (1998) show ev- 
idence for hot gas emerging from the inner regions to the 
north, which the authors interpret as evidence for a hot, 
expanding superwind fluid. 

The present Very Large Array (VLA) 21cm observations 
of the neutral hydrogen (H i) in Arp 299 were conducted 
as part of a larger program to study the closest lumi- 
nous IR galaxies in the IRAS Bright Galaxy Sample of 
Sanders et al. (1991). Arp 299 is both the most IR lumi- 
nous member of this sample within a distance of 50 Mpc 
{Lju = 8.1 X IO^^Lq), and the closest system with an IR 
luminosity greater than 5xlO^^L0. It is therefore an im- 
portant target to study for insight into the causes of such 
luminous starbursts. The purpose of the present study 
is to map any extended atomic hydrogen in order to con- 
strain the encounter geometry and the Hubble types of the 
progenitors. 

The paper is organized in the following manner. In ^ 
we briefly describe the observations. In §||we present the 
observations separately for the inner regions, tidal regions, 
and for the H i absorption and a new H i detected compan- 
ion. In §11 we use the observations to explore the probable 
Hubble types of the progenitors, the encounter geometry, 
the possibility of self-gravitating entities within the tail, 
the conditions for an IR luminous phase, and the expla- 
nation for the differing H i and optical tidal morphologies. 
We summarize our main conclusions in 

2. OBSERVATIONS 
2.1. HI Observations 

The VLA H i spectral line observations were obtained 
as part of a study on the outer gas dynamics of IR lumi- 
nous mergers (Hibbard & Yun 1996 and in preparation), 
and the data reduction techniques are fully described in 
that paper, except for the continuum subtraction which 
is described below. The details of these observations are 
tabulated in Table H^. Briefly, the data consist of a 3.5 
hour observation with the VLA in the C-array configura- 
tion (full-width at half maximum resolution 0pwHM^^5") 
and a 3 hour observation with the VLA in the D-array con- 
figuration (0_FVFffM~45"). The data were combined in the 
UV plane to form the C-f D array data that are used in the 
remainder of this paper. The correlator mode was chosen 
to provide 63 spectral channels with a channel spacing of 
10.5 km s~^ and total bandwidth of 660 km s~^ centered 
at 3080 km s^^. The data were reduced using standard re- 
duction procedures in the Astronomical Image Processing 
System (AIPS). 

The observations reveal broad H i absorption against the 
disk-wide starburst and central radio continuum sources 
(see also Dickey 1982; Baan & Haschick 1990 hereafter 



BH90). This absorption was observed in channels 5-55 
out of the 63 channel data cube (of which only channels 
2-58 are useful). As a result, we could not achieve an opti- 
mal continuum subtraction by using line free channels on 
either end of the bandpass, as is usually done (e.g. Rupen 
1998). Instead, we mapped and "cleaned" the continuum 
image and subtracted the brightest 17 "clean components" 
from the UV line data, which accounts for approximately 
half of the total continuum fiux of 820 mJy falling within 
the primary beam (FWHM=30'). The line data were then 
mapped and cleaned in AIPS, resulting in a three dimen- 
sional H I data "cube" , with Right Ascension and Decli- 
nation along the first two coordinates and velocity along 
the third axis. 

These data still contained many continuum sources, but 
were relatively free of the deep sidelobes from the brighter 
sources. The residual continuum was removed in the map 
plane using a spatially variable continuum baseline. The 
H I tidal features appeared over a narrow range of veloc- 
ity, and the residual continuum was removed from these 
regions by fitting a first order polynomial to a large num- 
ber of channels on either side of this range. The emission 
from the NW disk appears in channels 4-26, and the resid- 
ual continuum was subtracted from this region by fitting 
to the continuum in channels 34-54. The emission from 
the SW disk appears in channels 24-50, and the resid- 
ual continuum was subtracted from this region by fitting 
to the continuum in channels 4-20. Finally, the residual 
continuum was removed from the regions showing H i ab- 
sorption by fitting to the continuum in channels 2-4 and 
52-58. 

The data were mapped using different weighting func- 
tions to achieve varying resolutions and sensitivities. A 
high resolution data cube was made using a "Robust" 
parameter (i?; Briggs 1995) of R=0, which gives more 
weight to longer baselines and hence to smaller angular 
scales. This cube obtains a surface brightness sensitivity 
of 0.35 mJy beam~^ at a resolution of Ofwhm— ^7" xl5" . 
This corresponds to a column density limit (2.5ct) of 
4xl0^^cm~^ averaged over the beam width of 4.0x3.4 
kpc^. A more sensitive intermediate resolution data cube 
was made with R = 1, giving a resolution of 22"x20" 
(5.1 x4. 7 kpc^) and a column density limit of 2x lO^^cm^^. 
To further increase sensitivity to extended low column den- 
sity gas, a low resolution data cube was made by convolv- 
ing the R=l data cube to a resolution of 35" (8.1 kpc), 
reaching a detection limit of 1 x lO^^cm^^. These data will 
be referred to in the following as the high, intermediate, 
and low resolution data, respectively. 

2.2. Optical Observations 

The VLA observations reveal a gaseous stream of mate- 
rial extending to a radius of 9' that had no known optical 
counterpart on existing images. As such, this region was 
targeted for deeper optical observations with the Univer- 
sity of Hawai'i 88" telescope at Mauna Kea Observatory. 
The observational parameters are listed in Table |^. The 
Tek 2048 CCD was used with the f/10 re-imaging optics, 
giving a plate scale of 0.22" pixel"^ and a field of view of 
7'.5. Three overlapping 600 sec i?-band images of the re- 
gion of interest were obtained in January 1995. In January 
1997 we obtained a 600 sec i?-band image centered on the 
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main body and two overlapping 900 sec _B-band images 
of the H I tail. The total imaged region covers an area of 
8f7x 12.'4. The conditions were photometric on both dates, 
and the seeing was The data were calibrated via ob- 

servations of Landolt UBVRI standards (Landolt 1983) 
observed on the same nights, with zero point errors (Icr) 
of 0.01™ in B and 0.03™ in R. 

The images were flattened and combined using the tech- 
niques described in Hibbard & van Gorkom (1996, here- 
after HvG96), and the final mosaics are flat to better 
than one part in 500. The combined image was trans- 
formed to the World Coordinate System (and thereby to 
the same reference frame as the radio images) by register- 
ing to an image of the same area extracted from Version 
II of the Digitized Sky Survey^ (XDSS), obtained from 
the CD jukebox at the Canadian Astronomy Data Cen- 
tre.^ The registration was accomplished by referencing 
the location of ^20 stars in common on both images us- 
ing the "koords" program in the Karma software package 
(Gooch 1995) to perform a non-linear least squares fit for 
the transformation equations. The plate solution so found 
is better than a fraction of a pixel (0.22"), and the over- 
all registration should be as good as that of the northern 
portion of the Guide Star Catalog, which is estimated to 
be O'.'S (Taff et al. 1990). 

Finally, a deep optical image was constructed by bin- 
ning the pixels with the lowest light levels 9x9 to achieve 
a limiting surface brightness of /iB=28.2 mag arcsec"^ and 
^7?=27.1 mag arcsec"^ (2.5cr; Table ||). 



distribution in the inner regions. We have also labeled a 
minor axis clump of H i and a western plume for reference 
in the following discussion. 

To the right and below this panel are position-velocity 
maps constructed by summing the emission over declina- 
tion (lower panel) or right ascension (right panel). Fig. || 
presents the high resolution channel maps of the inner disk 
alone, after Hanning smoothing by a factor of two in ve- 
locity. In both of these maps, absorption is indicated by 
white dashed contours. 

The inner H i disk was previously mapped with the VLA 
in the C-array configuration by Stanford & Wood (1989). 
As discussed by these authors, the kinematic signature 
seen in Figs. || & || is representative of a rotating disk, 
with the NW half of the disk receding and the SW half 
approaching. The NW component of the disk reaches a 
radius of 25 kpc and contains 2.9x10^Mq of Hi. The SE 
component of the Hi disk reaches a radius of 21 kpc and 
contains a total of 2.8xlO^M0 of atomic hydrogen.® The 
disk spans the velocity range from 2880-3390 km s~^, with 
its kinematic center at 3135 km s~^. 

This disk was also mapped by the VLA in the D- 
configuration by Nordgren et al. (1997) using a coarser 
velocity resolution (by a factor of four) and a larger band- 
pass (by a factor of two) than the observations presented 
here. Their observations suggest that the emission from 
the NW disk continues to 3550 km s~^ (see their Fig. ^), 
and therefore it is likely that we have missed some disk 
emission due to our smaller velocity coverage. The OH 
maser emission from Arp 299 (Baan 1985; BH90) suggests 
that there may be nuclear gas at velocities as high as 3650 
km s^^. 

The H I disk appears to be physically associated with the 
optically distorted disk of IC 694. This view is strongly 
supported by both the CO and Ha kinematics ma pped at 
higher spatial resolution, as will be discussed in §4.1 be- 



3. RESULTS 

The luminosity for the main galaxies was calculated out 
to the /^B=26 mag arcsec"^ isophote, and the correspond- 
ing B magnitude for the Arp 299 system (tob=12.31) 
agrees within 0.01™ to that given by Mazzarella & Boroson 
(1993). These and other global parameters for the Arp 299 
system are given in Table ^ and the entries are described 
in the notes to that table. In the following, we will discuss 
the kinematics and morphology of the inner disk first, and 
then that of the outer tidal features. We also discuss the 
central H i absorption and report the discovery of a small 
dwarf companion to Arp 299. 

3.1. Inner Gas Disk 

The disk H i kinematics are illustrated in Figures || & 
1^. Fig. ^ presents several 2-dimensional representations of 
the high resolution H i data, with the emission summed 
over the third dimension. The upper left panel shows the 
integrated intensity map in solid contours (lowest contour 
corresponds to A''^f/'^5x lO^^cm"^) upon the i?-band im- 
age. This map emphasizes higher column density material 
at the expense of the diffuse lower column density H i. We 
have therefore provided a dotted contour to indicate the 
Nhi— 2x10^^ cm"^ contour of the intermediate resolution 
data, in order to give a more complete picture of the H i 

®The Second Palomar Observatory Sky Survey was made by the California Institute of Technology with funds from the National Science 
Foundation, the National Geographic Society, the Sloan Foundation, the Samuel Oschin Foundation, and the Eastman Kodak Corporation. 
The Oschin Schmidt Telescope is operated by the California Institute of Technology and Palomar Observatory. 

'^The Canadian Astronomy Data Center is operated by the Dominion Astrophysical Observatory for the National Research Council of 
Canada's Herzberg Institute of Astrophysics. 

^These masses do not include the contributions due to gas in the western plume or minor axis clump discussed later. 



low. There are, however, emission features which deviate 
from this simple disk motion. There are 4.6x10®Mq of 
gas in a western extension to the NW portion of the disk, 
and which appears kincmatically associated with that disk 
(i.e., appears in the same velocity channels in Fig. but 
which is spatially associated with a faint plume of starlight 
reaching to the SW (labeled "western plume" in Fig. ^, see 
also Fig. |l]). There are also 1.5x10®Mq of gas in a clump 
lying along the SW minor axis of the disk (labeled "minor 
axis clump" in Fig. |^). If this minor-axis material were as- 
sociated with the disk of IC 694, it should appear in only a 
few spectral channels. Since it appears over a broad range 
of velocities (from 3050—3200 km s'^ in Fig. |), we sug- 
gest instead that it is associated with NGC 3690 (see §4.1). 
There is also some emission associated with the base of the 
tidal tail (seen in the right panel of Fig. ^ towards the north 
and in the channel maps from 3092 — 3135km s~^), which 
is seen more clearly in the low resolution maps discussed in 
the next section. Finally, there are some H i features that 
project close to the compact spheroidal galaxy MCG+10- 
17-2a (Fig. H, channels at 3306 — 3371 km s~^). Since these 
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features appear at velocities that are higher than the ve- 
locity of MCG+10-17-2a itself (V^e/=3100 km s-\ FairaU 
1971, Sargent 1972), and close to the velocites of the Hi 
within the underlying disk of IC 694, it seems likely that 
they are due to disk material that was perturbed by a 
recent close passage of the spheroidal. 

3.2. Outer Tidal Morphology 

Our VLA observations revealed a previously unknown 
tidal tail reaching 124 kpc in radius to the north. This 
discovery was initially reported in Hibbard & Yun (1996). 
Subsequently, VLA D-array observations of Arp 299 from 
May 1995 were published by Nordgren et al. (1997), con- 
firming both the H i and optical tidal tail we report on in 
this paper. 

In Figure |^ we show the H i integrated intensity map 
from the intermediate resolution data, with position- 
velocity profiles plotted along either axis as in Fig. ^ 
The tail appears as two parallel filaments separated by 
about 20 kpc (labeled "inner filament" and "outer fila- 
ment" in Fig. H), extending northward and connecting 
onto or merging into a dense clump of gas (labeled "N 
clump" in Fig. 0). The inner filament originates near the 
minor axis of I(T 694 (see also the northernmost clump in 
Fig. 0), appears broken at the lowest contour drawn in 
Fig. ff (7V^f/~2.5x lO^^cm"^) after which there appears an 
Hi density peak {NHi,peak 2xl0^°cm~^, labeled "Hi 
knot" in Fig. ^), and the filament continues northward. 
This inner filament contains a total of 4.4xlO^M0 of neu- 
tral atomic hydrogen. The outer filament originates to the 
west of the first (see also the second northernmost clump 
in Fig. is more complete and of a higher average col- 
umn density than the inner filament (A^ff/~lxl02°cm-2), 
and contains twice as much Hi {Mhi=9x10^Mq). The 
filaments join at the N clump, which contains 2x lO^Af0 of 
H I. This clump has an irregular column density distribu- 
tion with peaks as high as 3xl0^°cm~^. The kinematics 
of these features will be discussed in the next section. 

The tail reaches a maximum projected radial separation 
of 8.'6 (124 kpc; measured from source A to the outer con- 
tour in Fig. and extends 12.'7 or 180 kpc as measured 
along its length (measured from source A along the peak 
H I column density contour of the outermost feature) . It 
contains a total of 3.3x10^ Mq of atomic hydrogen, which 
is at the upper end of the tidal gas content of other ma- 
jor disk-disk mergers (HvG96, Hibbard & Yun in prepa- 
ration). Since less than half of the progenitor disk mate- 
rial can be raised into a tail, and most of this falls back 
quickly onto the remnant (Hibbard & Mihos 1995), this 
amount of tidal H i requires a very H i rich progenitor 
{Mhi 10^°A/o). This point wiU be explored m more 



detail in §4.2 



The deep optical imaging reveals a faint stellar tail ex- 
tending in the same direction as the H i features (see Fig- 
ures H, H & 0). The optical tail appears to emerge in the 
north near the minor axis of the system and curve slightly 
to the NE to a radius of 9'. The width of the optical tail 
appears remarkably constant at '^l' (14 kpc), although 
there may be some widening in the i?-band image. The 
observed width of the tail can be understood as the result 
of the dispersal of outer disk stars with a typical velocity 
dispersion of 20 km s~^(e.g. van der Kruit 1988; Bottema 



1993) over an interaction age of 750 Myr (§4T). At very 
faint light levels the tail appears to loop south and west 
from the NE tip towards the bright star labelled "Star 2" in 
Fig. I^a, with perhaps a very faint extension directly south 
of this star. While scattered light from the stars labelled 
1 & 2 in this figure make the precise morphology of the 
faint optical light uncertain, these features have the same 
morphology in both the B- and i?-band greyscale images 
of Fig. H, while the scattered light properties are differ- 
ent between the bands (see e.g. the difference between the 
faint scattered light associated with Star 3 in Figs. ||a & 
|5|3), and we feel confident that the faint optical features 
are not artifacts. 

Upon first examination the optical and H i tidal fea- 
tures appear similar. But when the data are overlaid, as in 
Figs. II& the differences in the distribution are striking: 
the inner H i filament roughly aligns with the optical tail, 
while the outer parallel H i filament has no correspond- 
ing optical feature down to the faintest levels measured 
mag arcsec"^, ^ji^27.5 mag arcsec"^). 

Even within the optical tail, there is a poor correlation 
between optical and gaseous density peaks. If anything, 
the gas peaks seem to be displaced a bit to the west of 
the optical peaks. Neither the gap nor the knot in the 
inner H I filament appear to coincide with any discrete op- 
tical feature, although the optical surface brightness drops 
by ^1 mag arcsec^^ just north of the knot (Fig. |^). The 
displacements are even more extreme in the northernmost 
regions of the tail, where the gas and starlight appear to be 
anti-correlated: to either side of the optical tail there are 
regions of high gas column density (Af/f/'^3x lO^^cm"^ to 
the west, iVij/~2xlO^"cm^^ to the east), whereas the gas 
within the tail has iVjj/~lxlO^°cm~^. At its northern- 
most extent, the optical tail has no associated H i, and ac- 
tually appears to "hook" towards the east, exactly around 
the northern H i contours in Fig. ^a. 

This anti-correlation suggests the possibility that dust 
associated with the H i gas is playing a role in shaping the 
optical tidal morphology. However the required amount of 
extinction (>1 mag at Nhi— 2xl0^°cm~^) is an order of 
magnitude larger than measured for Galactic H i (Bohlin 
et al. 1978). Since tails are drawn from the outer, presum- 
ably dust-poor regions of disks, such a high gas-to-dust 
ratio seems unlikely (see also Alton et al. 1998). Evidence 
for just such a high gas-to-dust content in tidal gas may be 
indicated by observations of red globular clusters behind 
an HI tidal stream in NGC 5018 (Hilker & Kissler-Patig, 
1996). Still, we feel that this explanation is ruled out for 
Arp 299 by the fact that the optical tail is actually bluer 
north of the Hi knot, where the Hi column density is 
highest, compared to the starlight in the "gap" in the in- 
ner Hi filament, where the Hi column density is lowest 
(Fig. ^). One would expect the opposite trend if there 
was a high gas-to-dust ratio in this tidal gas (however see 
Witt, Thronson & Capuano 1992). The B - R colors of 
the stellar tail range from 0.3-1.8, similar to range found 
in the outer stellar disk of late type galaxies (de Jong 1995, 
Ch. 4) and in other tailed mergers (Schombert et al. 1990; 
Hibbard et al. 1994; Hibbard 1995). 

The total stellar content of the optical tail was es- 
timated by masking the optical images below ^b= 28 
mag arcsec~^and iir= 27 mag arcsec"^ and replacing re- 
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gions con tain ing; point sources (presumably background 
sources, § |4.5| ) and the optical halos of Stars 1& 2 in 
Fig. ||a with adjacent background values. The total op- 
tical light in the northern tail is Lb = 1.6 x IO^Lq b and 
Ln = 1.2 X IO^Lq^h, corresponding to 4% and 2% of the 
total B- and i?-band light of the system, respectively^ (Ta- 
ble H) . These percentages are on the lower end of the values 
found in peculiar systems (Schombert et at 1990; Hibbard 
et al. 1994; HvG96). 

The total H i mass-to-blue light ratios (in solar units) 
for all the northern tidal features is Mhi / Lb= 1-8. There 
is quite a variation within the tail, with Mhi /Lb "^0-1 
along the optical tail, and ranging from 2 to over 10 for 
the outer filament and eastern regions of the N clump. The 
faint plume of light extending 45 kpc in radius to the west 
of Arp 299 contains = 1 x 10^ Lq^b- Associating the 
western Hi extension to this feature gives Mhi / Lb=0.5. 

3.3. Outer Tidal Kinematics 

The kinematics of the tidal regions are illustrated by 
the integrated position- velocity profiles in Fig. ^, intensity 
weighted velocities and velocity dispersions in Fig. @c&d, 
individual channel maps in Fig. and two position- 
velocity slices through either H i filament in Fig. ^. Both 
Fig. ^ and Fig. ^ are constructed from the low resolution 
data cube, and the lowest plotted contour corresponds to 
a column density of IxlO^^cm"^. The channel maps are 
only plotted over the velocity range containing tail emis- 
sion. 

The kinematics immediately illustrate two important 
points. The first is that the tidal kinematics are smooth, 
continuous and single valued, as expected for tidal features 
(e.g. Barnes 1988; Hibbard et al. 1994; Hibbard & Mihos 
1995). The second is that the two H i filaments s hare the 
same kinematics. This is best seen in Fig. ||, in which the 
H I emission associated with parallel regions along the fila- 
ments appears in identical channels all along the tail. This 
rules out the possibility that the two filaments are separate 
tidal tails that are simply projected near each other; they 
must be different regions of the same kinematic structure. 
The figures also show that the N clump is continuous with 
the filaments and not a distinct entity. Therefore, the N 
clump and two filaments form a single contiguous feature 
in both space and velocity. 

A more detailed comparison of the filament kinematics 
is shown in Fig. ^, where we have plotted position- velocity 
profiles along both. This figure shows that the base of the 
outer filament has a slightly higher mean velocity than 
the base of the inner filament, which may help in decid- 
ing the relative morphologies, although we await detailed 
kinematic modeling before deciding how to interpret this 
information. Fig. || illustrates a few more interesting prop- 
erties of the filaments. In Fig. ^b we see that the apparent 
"gap" in the integrated intensity map of Fig. ^ is not de- 
void of H I: weak emission which did not pass the threshold 
used when making the moment maps bridges this gap (see 
also the panel at 3150 km s^^ in Fig. We also see that 
the "knot" within the inner filament has an anomalously 
large velocity width (see also line- width map in Fig. |^): 
the majority of the tail has a velocity dispersion of 7-9 

^These numbers increase to Lb = 2.1 X W^Lq^b (5%) and Ln = 
mag arcsec"^ lower surface brightness threshold is used. 



km s^^, which is typical of tidal tails in general (Hibbard 
et al. 1994, HvG96), while the knot has a dispersion of 
over twice this (20 km s^^). Fig. |^ shows that there are 
two additional velocity dispersion peaks along the south- 
ern edge of the N clump, with dispersions of 17 km 
and 13.5 km compared to an average dispersion of 10.5 
km s^^ within the N clump. None of these regions corre- 
sponds to any noticeable features in the optical images or 
color maps, which argues against their being related to the 
putative tidal d war f galaxies, a question we investigate in 
more detail in §4.5. 



Finally, there is a high velocity envelope of material just 
north of the knot (indicated in Fig. ||b). A similar enve- 
lope of material is seen at parallel locations in the outer 
filament (Fig. ^). The channel maps show that this gas 
lies exactly along the faint optical light that stretches from 
the end of the tail towards Star 2 in Fig. ||a (see panels at 
3171-3161 km s^^ in Fig. ^), while the majority of the fil- 
ament gas has lower velocities and lies north of this light 
(see panels at 3098-3140 km s^^ in Fig. ^). These are 
probably important clues to the origin of these filaments, 
although we are at a loss to say what. 

3.4. Central Absorption 

High resolution (~5") VLA radio continuum observa- 
tions show that half of the total the 1.4 GHz radio con- 
tinuum emission lies within 5" of regions A, B and C-C 
in Fig. |l|, with the other half of the flux distributed be- 
tween these regions and throughout the disk of IC 694 
(Gehrz, Sramek & Weedman 1983). Even higher resolu- 
tion {9fwhm=^''7) VLA A-array observations of the Hi 
absorption against these continuum sources are presented 
in BH90, along with observations of OH megamaser emis- 
sion. These observations show that H i absorption is seen 
against both the compact continuum sources as well as 
against the extended continuum features. Our continuum 
profiles are quite similar to those published by BH90, with 
subtle differences that are likely due to our larger beam, 
which allows surrounding H i emission to fill in some of the 
absorption. We will therefore use the results of the higher 
resolution H i absorption observations in concert with our 
H I emission observations to set some rough limits on the 
spin temperature (Tgpin) of the H i gas. Along the way, we 
will estimate the absorbing column of gas in front of the 
continuum sources, from which we will derive an approxi- 
mate correction to the total H i mass. 

BH90 found that the H i absorption is dominated by 
a broad component {AV ^120 km s^^) at V';ie;=3160 
km s~^, and a narrow component {AV ~70 km s~^) at 
Vhe/=3050 km s~^. Due to the near correspondence of 
these velocities with the systemic velocities of IC 694 and 
NGC 3690 (3080-3180 km s'^and 3000-3080 km s'^, re- 
spectively; BH90 and reference therein) and because these 
two components are found in the absorption spectrum 
against all of the radio continuum sources, BH90 deduce 
that they are due to outer disk material from both IC 694 
and NGC 3690 lying in the foreground rather than near- 
nuclear gas clouds. This allows us to interpolate our H i 
column density maps across the absorbing regions. This 
column density can then be used together with the ratio of 

1.9 X lO^-LQ.ij (4%) if the background point sources arc left in and a 0.5 
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the neutral gas absorbing column density to the spin tem- 
perature (Nhi /Tspin) given by BH90 to constrain Tgpin- 

From the column densities given in Fig. 2, we estimate a 
foreground column density of N^i «l-2x lO^^cm"^ due to 
the disk of IC 694. For the corresponding absorption com- 
ponent (i.e. the broad component at V/ie;=3050 km s^^), 
BH90 find Nni/T^pin -S-SxlO^^cm-^ fg K-\ where fg 
is the fraction of the radio continuum source covered by 
absorbing gas. Since BH90 find this absorption feature 
in the spectra of each continuum source in Arp 299, fg 
must be near unity, yielding a spin temperature of 12-70 
K. This is similar to the range found in our galaxy (34-74 
K; Kalberla et al. 1985). The lower range of temperatures 
are similar to the very low values of Tspin (Tgpin < 50 
K) found in the Magellanic stream (Kobulnicky & Dickey 
1999) and the nearby galaxies NGC 247 and NGC 253 
(Dickey, Brinks & Puche 1992). 

Using this range of spin temperatures and the mea- 
sured N HI /Tspin values from BH90 for the low veloc- 
ity absorption component (V;ie!=3050 km s~^, Nni/Tspin 
= l-2xl0^^cm~^ K~^) we deduce a column density of 
1020 _ io2icm-2 due to gas from NGC 3690. Adding these 
two contributions^" gives an expected column density of 1- 
3xl0^^cm^'^ in front of the absorption region. Doubling 
these values to account for emission behind the continuum 
sources we conservatively estimate an average column den- 
sity in the absorbing region of 2xl0^^cm~^, corresponding 
to a "missing" Hi mass of 2.2x10®Mq. This value is used 
to arrive at the "absorption corrected" Mhi value tabu- 
lated in Table ||. This total is very approximate, and could 
easily be 2-3 times as large (see §4.2). 



3.5. HI Detected Companion 

The H I observations uncovered a previously unknown 
and uncataloged companion to Arp 299 containing 8 x 
10'' M0 of H I. This companion lies ~10' west and 3' south 
of Arp 299 (projected separation of r=140 kpc) and has 
a systemic velocity of -1-145 km s^^ with respect to Arp 
299. This region was not covered in our optical observa- 
tions. An amorphous optical counterpart is found at this 
location on the XDSS image. The Hi column density is 
contoured upon the XDSS image in Figure |o[ and channel 



maps are shown in Figure 11 



We attempt to derive an optical luminosity for this 
companion from the POSS-E plate of Vers ion I of the 
DSS, using the photometric solution given at http://www\ 
gsss.stsci.edu/dss/photometry/poss_e.htTn\ . However, the 
number of counts from this source (log [counts] = 4.52) 
is lower than the levels over which calibration has been 
reliably determined (log[counts] = 4.8, mposs-E < 16 
mag). Extrapolating the faint end of the above curve to 
the observed number of counts gives a somewhat uncertain 
estimated magnitude of mpQgs-E^^^ mag. Assuming an 
average value oi B — R — +1.0 (e.g. van Zee, Haynes, & 
Salzer 1997), the absolute B magnitude of the Hi com- 
panion is about Mb = —15.4. 

The resulting global properties for the companion are 
given in Table 0. The optical appearance and H i proper- 
ties of this system are similar to those of dwarf irregulars 

^"There are additional absorption components seen in the H I absorption spectra obtained by BH90, but their association with OH omission 
features at similar velocities places them very close to the nucleus of IC 694 with a small covering factor, and they will therefore contribute 
very little to the total H I mass. 



(Hoffman et al. 1996), and it is likely to have been a dwarf 
irregular companion of one of the progenitors. For compar- 
ison, it is about 20 times fainter and 30 times less gas-rich 
than the LMC (TuUy 1998). 



4. DISCUSSION 

Although the differences between the optical and 
gaseous tidal morphology are bizarre, the tail kinemat- 
ics and general morphology are very similar to the 
gas-rich tails of more classically understood mergers 
(e.g. Schweitzer 1978; van der Hulst 1979; Hibbard et al. 
1994; HvG96), and the two main bodies certainly appear 
to be separate disk systems that have been strongly per- 
turbed. It therefore seems that the tidal hypothesis — 
that these features arise as a result of two spiral galaxies 
falling together under their mutual gravitational attrac- 
tion and merging due to the effects of dynamical friction 
(e.g. Toomre & Toomre 1972 hereafter TT72; Barnes & 
Hernquist 1996) — still provides the simplest explanation 
for the many peculiarities in this system. 

Given this premise, we now aim to understand the gross 
characteristics of the encounter: what was the general en- 
counter geometry (j j4.l| ), and what we re the approximate 
Hubble types of the progenitors (§4.2)? We explore pos- 
sible explanations for the differences between the gaseous 
and stellar tidal features in j |4.3| , although much of this 
discussion is deferred to a separate paper (Hibbard, Vacca 
& Yun 1999, hereafter HVY99). In §4.4 we examine what, 
if anything, the state of this system tells us about the con- 
ditions under which super starbursts are initiat ed, leading 
to an ultraluminous infrared phase. Finally, in §4.5 we dis- 
cuss the lack of putative tidal dwarf galaxies in the outer 
tidal features. 



4.1. The Nature of the Encounter 

The tail kinematics do not themselves determine the 
spin geometry of their progenitor. This is because the line- 
of-sight velocity fields of tidal features are often dominated 
by projection effects. This is particularly true of long fil- 
amentary tails, for which the space velocities will be ori- 
ented primarily perpendicular to our line-of-sight (leading 
directly to their large apparent lengths). Instead, we use 
the morphology of the tail to deduce the relative spin ge- 
ometry between its progenitor disk and the orbital plane. 
In particular, only disks with a prograde spin geometry 
(i.e., a spin angular momentum roughly aligned with the 
orbital angular momentum of the two merging galaxies) 
can raise such long, drawn out features (TT72, Barnes 
1988), and the progenitor of the northern tail must have 
had just such a prograde spin orientation. 

Further, since tidal tails remain close to the original spin 
plane of the disk from which they are extracted (TT72), 
and the observed tail in Arp 299 emerges at nearly right 
angles to the disk of IC 694, we conclude that IC 694 could 
not have given rise to this feature. The tail must therefore 
have arisen from NGC 3690. Similarly, the lack of any 
narrow drawn-out features near the plane of the disk of IC 
694 suggest that this system experienced the encounter in 
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a very highly inchned or retrograde sense. Such encounters 
provoke a much milder kinematic response in a disk (e.g. 
TT72; White 1979; Noguchi 1991), accounting for the sur- 
vival of the disk morphology and kinematics of IC 694 so 
late into the encounter. This spin geometry agrees with 
the kinematics of the inner regions derived from Ha emis- 
sion line kinematics, which show that IC 694 and NGC 
3690 rotate in opposite directions (Augarde & Lequeux 
1985; Hibbard, Bland-Hawthorn & TuUy, in preparation). 
Further, the systemic vel ociti es for the two nuclei (Vhei ~ 
3135 km s~i for IC 694 [§|I1 and Vhei « 3040 km s'^ for 
NGC 3690 [§ ^.4[ ) give them the proper sense of motion 
for the inferred spin geometry — the nuclei rotate about 
each other opposite the direction that the disk of IC 694 
rotates. 

The morphology and kinematics of the ionized gas in 
Arp 299 has been mapped with Fabry-Perot interferome- 
try (Hibbard, Bland-Hawthorn & TuUy, in preparation), 
showing that there is spatial and kinematic continuity be- 
tween the outer H i disk all the way inward to source A. 
This implies that the entire disk structure belongs to IC 
694; i.e. the H i is not distributed in a settled disk in rota- 
tion about both systems, as suggested by Stanford & Wood 
(1989). Further, these observations reveal that NGC 3690 
and sources B &c C-C are kinematically distinct from 
the IC 694 disk, as expected from the inferred prograde- 
retrograde spin geometry. 

From the redshifted velocities along most of the outer 
tail, and the fact that the tail should be predominantly 
in expansion along most of its length (Hibbard & Mihos 
1995), we deduce that the tail curves away from us and 
should connect to the main bodies from behind. In this 
case, the decrease in velocities seen at the very base of 
the tail in Fig. ^3 (3rd arrow from the bottom in both 
figures; see also northern region of upper right panel in 
Fig. indicates that the material at the base of the tail 
has reached its apocenter and is falling back to smaller 
radii (e.g., Hibbard et al. 1994; Hibbard & Mihos 1995). 
Similarly, the broad range of velocities seen for the minor 
axis clump of H i (Fig. ^ suggests that this material is 
part of the northern tail connecting back to its progenitor, 
NGC 3690. If there is more Hi associated with the disk 
of NGC 3690, it would be difficult to disentangle from the 
emission due to the SW portion of the disk of IC 694 and 
the central HI absorption. 

Finally, the western plume may be disk material that 
was perturbed during an earlier pericentric passage of the 
two galaxies about each other or by the passage of the com- 
pact NW spheroidal (MCG-)-10-17-2a). Since tidal mate- 
rial tends to lie within the plane it occupies at the time 
it is raised (TT72), the distribution of this gas should re- 
flect the plane of the disk at an earlier stage of the merger, 
perhaps at an early pericentric approach. The fact that 
this plane is different from the present plane suggests that 
the angular momentum of the disk was changed during 
the encounter, evidence for which is also seen in the Ha 
Fabry-Perot maps at smaller radii. 

These observations show the power of sensitive H i and 
deep optical imaging of peculiar systems. It has long been 
believed that Arp 299 is in a very early stage of its inter- 
action, given the relative youth of the on-going starburst 
(^20 Myr; Vacca et al. in preparation; Meurer et al. 1995; 



Stanford 1989; Nakagawa et al. 1989; Augarde & Lequeux 
1985). The apparent lack of tidal features in their H i ob- 
servations led Standford & Wood (1989) to conclude that 
the outer tidal H i had quickly settled into a disk config- 
uration early in the encounter. The discovery of the 180 
kpc tidal tail obviously changes this conclusion. It also 
requires the interaction to be quite old. Taking a typical 
rotation speed of 240 km s~^ from the inner disk kinemat- 
ics and saying that this is the maximum velocity available 
to the tidal debris, we calculate that the tail has taken at 
least 750 Myr to form. Therefore the dynamical age of the 
merger is substantially larger than the age of the present 
starburst. Arp 299 is not unique in this sense (e.g. Mihos 
& Bothun 1998), and this cautions against using the age 
of the starbursts to age-date an interacting system, as is 
frequently done for the peculiar star-bursting objects at 
high redshift. Further, it is possible that Arp 299 experi- 
enced several short bursts of elevated star formation since 
the tail was first launched. In this case, the accumulated 
burst populations would not be co-eval, but would show 
a spread of ages, and the expected age signatures in the 
resulting remnant will be much less pronounced then is 
often presumed {e.g. Silva & Bothun 1998). 

4.2. The Nature of the Progenitors 

Now that we can use the H i kinematics to approxi- 
mately associate the different kinematic features with their 
progenitors, we attempt to divide the gas content and B- 
band luminosities between the two systems. We will then 
use these global quantities in concert with the statistical 
properties of normal Hubble types in an attempt to de- 
duce the approximate Hubble types of the progenitors. It 
should be noted that not only are there significant uncer- 
tainties inherent in apportioning the gas and light between 
the two systems, but there is also a notoriously large scat- 
ter in global properties along the Hubble sequence, and a 
further uncertainty due to conversion between gas phases 
and from gas into stars due to the on-going starburst. 
However the following should provide a very rough guess 
for the likely progenitors. 

As mentioned above, we associate the H i from 
the tidal tail and minor axis clump with NGC 3690 
(M///=3.5x lO^M©), and the inner disk and W plume of 
Hi with IC 694 {Mhi=6.2x10^Mq). Simple dynamics 
suggest that less than one-half of the disk can be ejected 
into a tidal tail, and that most of this quickly falls back 
onto the merging system (Hibbard & Mihos 1995). This 
suggests a total Hi content of 2-3 xMHi,tidai = (7- 
10) X lO^Afg for NGC 3690. However, this value is already 
more than the combination of the tidal H i and the absorp- 
tion correction (5.7xlO^M0). This discrepancy suggests 
that o ur a bsorption correction was probably too conserva- 
tive (§3.4). Rather than trying to account for gas that we 
have nonrm measure of, we simply associate all of the ab- 
sorption corrected gas (2.2x10^Mq) with NGC 3690 and 
remark that the H i masses for both IC 694 and NGC 
3690 are likely significantly underestimated. Using the 
12C0(1 — 0) observations of Casoli et al. (1999), we asso- 
ciate the "widespread" molecular gas and that near source 
A with IC 694, the molecular gas near source B with NGC 
3690, and split the molecular gas from the overlap region 
(C-C) between the two systems. Total gas masses are 



8 



Tidal Tail in Arp 299 



calculated by multiplying the Hi mass by 1.34 to account 
for He and adding this to the molecular gas mass. Finally, 
the individual i?-band luminosities are estimated by sum- 
ming the light over irregular polygons drawn around each 
system, and are therefore also very approximate. The re- 
sults of this division are given in Table |, and we use these 
values in concert with the following statistical properties 
of normal Hubble types: 

• Mh^/Mhi increases towards earlier Hubble types 
(Young & Knezek 1989 Young & ScoviUe 1991). 
The high values in Table |j suggests progenitor types 
later than Sab and earlier than Sc. They also argue 
against low surface brightness progenitors (de Blok 
1995) 

• Mhi/Lb increases towards later Hubble types 
(Roberts & Haynes 1994; GiovaneUi & Haynes 1990). 
The high values in Table ^ suggest types later than 
Sa (Wardle & Knapp 1986, Bregman et al. 1992), 
and later than Sb (Roberts & Haynes 1994). 

• The large mass of H i suggests types later than Sa 
(Huchtmeier & Richter 1988; Bregman et al. 1992) 
and earlier than Sc (Huchtmeier & Richter 1988; 
GiovaneUi & Haynes 1990). 

• The large values of A/^a and M gas /Lb suggest types 
later than Sa (Hogg et al. 1993). 

• The large values of MH2/LB favor types later than 
Sab (Lees et al. 1991, Bregman et al. 1992). 

These consideration suggest progenitor types of Sab-Sc 
for both systems - an earlier type for IC 694 and a later 
type for the NGC 3690. This view is supported by high- 
resolution (1") NIR imaging of the brightest sources at 
2.2 /im by Wynn-Wilhams et al. (1991). In that study, 
the light profiles of the brightest NIR sources within Arp 
299 were compared with stellar profiles taken immediately 
before and after the on-source exposures. It was found 
that source A is clearly resolved, with a FWHM of 460 pc 
(correcting to our adopted distance of 48 Mpc), while the 
point sources within NGC 3690 are found to have cores 
that are indistinguishable from stars, implying that they 
have diameters of less than 230 pc at 2.2/im. These con- 
siderations suggest a more extended bulge for IC 694 than 
for NGC 3690, lending further support to an earlier type 
for the former and a later type for the latter. We therefore 
suggest progenitor types of Sab-Sb for IC 694 and Sbc-Sc 
for NGC 3690. 

In summary, we believe that Arp 299 results from a 
prograde-retrograde encounter between a gas-rich Sab-Sb 
(IC 694) and Sbc-Sc (NGC 3690) that occurred 750 Myr 
ago. We estimate that the system will fully merge on a 
timescale somewhere between a crossing time (d/vmax) 
and an orbital period {'nd/vmax) at the present separa- 
tion d of 4.7 kpc, or sometime between 20-60 Myr from 
now (for a rotational speed of 240 km s^^). 

4.3. Explanations for Outer Tidal Morphology 



The differences between the gaseous and stellar tidal 
morphologies were very unexpected. Displacements be- 
tween optical and gaseous tails or bridges are not unheard 
of (e.g. NGC 4725/47, Wevers et al. 1984; Arp 295, HvG96; 
NGC 7714/5, Smith et al. 1997; NGC 2782, Smith 1994); 
nor are gaseous tidal tails or streamers with little, if any, 
starlight (e.g. NGC 2444/5, Appleton et al. 1987; M81, 
Yun et al. 1994); nor are stellar tails or plumes with no 
Hi (e.g. NGC 3921, HvG96; Arp 105, Due et al. 1997; 
NGC 2782, Smith 1994). There are even examples of much 
less extreme "bifurcated" gaseous tails, in which a single 
tidal tail appears to split near its base into parallel Hi 
filaments, both of which have similar kinematics (e.g. M81 
Yun et al. 1994; NGC 3921 HvG96; NGC 2535/6 Kaufman 
et al. 1997; NGC 4038 Hibbard, van der Hulst & Barnes, in 
preparation). What is unique in Arp 299 is that all these 
effects occur at once — there is a tail which is bifurcated 
along almost its entire length, with one filament gas rich 
and the other gas poor, and a striking anti-correlation be- 
tween the gaseous and stellar filaments where the filaments 
join. In this subsection, we investigate possible explana- 
tions for each of these characteristics in turn. 

4.3.1. Gas-Rich Outer Filament 

The high Mhi/Lb ratios found for the outer filament 
are rather extreme, with peaks >10 Mq Lq and a mean of 
~4 Mq Lq^. However they are similar to the values found 
in other gaseous tidal extensions (e.g. Arp 143 Appleton 
et al. 1987; NGC 2782 Smith 1994; NGC 7252 Hibbard 
et al. 1994). In these systems, the gas-rich extensions ap- 
pear to be due to the fact that Mhi/Lb is an increasing 
function of radius in spiral galaxies (Wevers et al. 1986)^^. 
Since the outermost (gas-rich) regions of the disk become 
the outermost regions of the resulting tidal tails (TT72), 
gas rich extensions to optical tails are to be expected. The 
higher dispersion of outer disk stars with respect to the gas 
(van der Kruit 1988) and the finite lifetimes of the most lu- 
minous stars will exacerbate this effect, leading to further 
increases in Mhi/Lb in the outer regions (HvG96). 

However, Arp 299 is quite different from the above men- 
tioned systems in that it is not just the outermost regions 
that are gas rich — the gas-rich filament appears to ex- 
tend all the way back to the main body, running parallel to 
the optical tail. That is, it appears to cover a very broad 
range of radii. So the first question we wish to address is 
whether such a feature can arise due to combination of dy- 
namical and projection effects which cause gas-rich outer 
regions to be projected adjacent to optically brighter inner 
regions. 

We believe that answer to this question lies in the dy- 
namical development of tidal tails, as demonstrated so 
vividly by TT72. These authors emphasized that tails 
are not linear structures, but actually two dimensional 
"ribbons" twisting through space. In particular. Fig. 2 
of TT72 shows that the outermost edges of the ribbon 
come from the outermost radii of the original disk, while 
middle regions come from intermediate radii, etc.; hence, 
the outer disk material forms a "sheath" around the inner 
disk material along its entire length. Therefore it is pos- 
sible that a lateral twist may cause the outer edge of the 

^^This is usually expressed as the well-known fact that H I disks often extend further than optical disks e.g. Bosma 1981, Wevers et al. 1986, 
Cayatte et al. 1994, Broeils & van Woerden 1994. 
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ribbon to lie in a different plane from the inner regions, 
i.e. for the most gas-rich regions to appear in a different 
plane from the less gas-rich, optically brighter regions. 

The models also show that the material from different 
radii form continuous structures in space and velocity, and 
hence that the resulting tidal features will necessarily be 
continuous in space and velocity as well (Barnes 1988, Hi- 
bbard & Mihos 1995). But while the parallel filaments 
of Arp 299 appear to be kinematically continuous, they 
appear morphologically distinct, or bifurcated. We sug- 
gest that the bifurcation arises due to a pre-existing warp 
in the outer regions of one of the disk. This suggestion 
comes from our recent attempts to match the tidal mor- 
phology and kinematics of the northern tail using N-body 
simulations (Hibbard et al. in preparation). In our prelim- 
inary trials (using only 4096 particles per galaxy), we can 
rather easily match either the morphology and kinematics 
of the outer tidal filament and N clump, or the kinematics 
of the inner filament and morphology of the stellar tail. 
However we have been unable to match both features si- 
multaneously with a single disk orientation and viewing 
angle. This leads us to suspect that the disk of NGC 3690 
was initially warped (as is often the case for disk galaxies; 
Bosma 1991). In this case, the outer, optically fainter re- 
gions will be pulled off along one plane (resulting in the 
outer tidal filament), while the regions from smaller radii 
will move primarily in a different plane (resulting in the 
inner filament). A similar explanation may account for 
the less extreme bifurcated structures seen in other tidal 
tails (e.g. M81 Yun et al. 1994; NGC 3921 HvG96; NGC 
2535/6 Kaufman et al. 1997; NGC 4038 Hibbard, van der 
Hulst & Barnes, in preparation). 

However, this projection effect does not explain the low 
H I content of the optical tail in Arp 299. Although the 
less extreme radii of spirals have a lower Mhi /L b than the 
outer regions, they still have significant quantities of H i. 
We therefore seek a separate explanation for this feature. 

4.3.2. Gas-Poor Inner Filament 

For the original disk to simultaneously give rise to a gas- 
rich outer filament and a gas-poor optical filament would 
require a special H i radial distribution — one in which the 
gaseous disk had a low column density within the stellar 
disk, but was very gas-rich beyond this. Such distribu- 
tions are indeed seen in some early type galaxies, espe- 
cially SBO's (van Driel & van Woerden 1991). However, 
two facts lead us to suspect that this is not a viable expla- 
nation for Arp 299. The first is that such systems usually 
have an order of magnitude less H i than observed in the 
tidal tail of Arp 299 (which itself is less than half of the 
original gas content). The second is that the starlight in 
such systems is dominated by an extended central bulge, 
which the NIR observations appear to rule out for NGC 
3690 (§0). 

Another factor that leads us to consider other explana- 
tions is that Arp 299 is one of a number of peculiar star- 
burst galaxies exhibiting curious displacements between 
tidal Hi and starlight. Further, the five systems which 
exhibit the most dramatic displacements (M82, Yun et al. 
1993; NGC 4631/4656, Weliachew et al. 1978; NGC 520, 
HvG96; Arp 220 Yun, Hibbard & Scoville in preparation, 
HVY99; and Arp 299) also host massive nuclear starbursts 



with associated powerful outflows that extend many kpc 
from the nuclear regions ( "superwinds" , e.g. Heckman, 
Lehnert & Armus 1993; Lehnert & Heckman 1996), and 
in each case the tidal H i shows a displacement or anti- 
correlation along the direction of the expanding hot super- 
wind fluid. This leads us to suspect a possible interaction 
between the expanding wind fluid and the tidal debris (see 
also Chevalier & Clegg 1985; Yun et al. 1993; Heckman, 
Lehnert & Armus 1993; Dahlem et al. 1996). 

This issue is covered in more detail in HVY99, and here 
we just discuss the specific details of Arp 299. In this sys- 
tem the evidence for an expanding superwind comes from 
the morphology and kinematics of the optical emission line 
gas (Armus, Heckman & Miley 1990) and the X-ray ob- 
servations of Heckman et al. (1998). In the latter, the soft 
X-ray emission from hot gas is oriented perpendicular to 
the harder X-ray emission from discrete x-ray sources as- 
sociated with the star burst. The hot gas emission is elon- 
gated roughly along the minor axis of IC 694, extending to 
a radius of order 25 kpc (Heckman et al. 1998). Since the 
observed X-ray emission represents just the hottest, dens- 
est portions of the adiabatically expanding wind (Wang 
1995), the full extent of the wind may extend well beyond 
this. 

The mechanical effect of such a wind may be evidenced 
in the increased column density and velocity dispersion of 
the Hi knot in the inner filament (r « 60 kpc), which 
might be due to the winds ram-pressure effects. The wind 
might also account for the high-velocity envelope of gas 
seen in Fig. |^ and the high velocity dispersion regions seen 
on the southern end of the N clump in Fig. However, 
we do not believe that the winds ram pressure has pushed 
the gas fully out of the inner filament and into the outer 
filament. Such action should significantly change the kine- 
matics of the gas, whereas the kinematics of the gas within 
both filaments is nearly identical. Instead, we suspect that 
the superwind affects the ionization state of the gas in one 
of two ways: either directly by shock ionization, or indi- 
rectly by clearing a sight-line from the tidal tail to the nu- 
clear starburst, allowing UV photons from the young stars 
to escape these regions and reach parts of the tail. Calcu- 
lations carried out in HVY99 suggests that both scenarios 
are in principle capable of affecting gas at the observed col- 
umn densities and distances, although there is insufficient 
evidence to discriminate between them. 

Whether the "missing" H i has been ionized or not 
is subject to observational investigation. The ex- 
pected emission measure in convienient units is EM = 
0.42cm-6pc(2^^|^)2(i2lE£)2, where Nhi is the 
missing column density of H i, and dL is the line-of-sight 
depth of the H i column. If the gas has a clumpier distribu- 
tion, then there should be some higher density peaks with 
emission measures well above these levels. These emmis- 
sion measures are within the capabilities of modern CCD 
detectors (e.g. Donahue et al. 1995 obtained 0.3 cm~^ 
pc in NGC 4631), and such emission would be well worth 
looking for. 

4.3.3. Gas/Light Anti- correlation in the N Glump 

The main problem with the above scenarios is that nei- 
ther can explain the H i and optical anti-correlation at 
the end of Arp 299 tail, whereby there is a local minima 
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in the H i column density where the optical tail crosses 
the N Clump, and local maxima to either side of this 
(Fig. ^). The expanding galactic wind or photoioniza- 
tion cone should have no idea where the optical tail lies. 
In HVY99, we suggest that the outer H i from the warped 
disk is already in a highly ionized phase, due to its low 
density and the intergalactic UV field. In this environ- 
ment, we suggest that the additional ionizing radiation 
due to evolved sources within the stellar tail itself, such as 
late B stars and/or white dwarfs, increases the ionization 
fraction of the local gas over that of the tidal gas without 
accompanying starlight. The lower neutral fraction is then 
seen as a drop in the neutral gas column density. A simple 
calculation is carried out in HVY99 that suggests that this 
process is in principle feasible, if the gas in the N clump 
has a line-of-sight thickness ^25 kpc. In thise case, the 

expected emission measure (^0.02 cm~^ pc) is well below 
what is observable with current technology. 



4.3.4. Outer morphology: conclusions 

These explanations are somewhat unsatisfactory since 
we need three separate mechanisms to explain the three 
puzzling characteristics: a warped H i disk to explain the 
parallel tidal filaments; a star burst wind or unusual initial 
H I distribution to explain the lack of H i within the op- 
tical filament; and an increased ionization fraction in the 
gas due to the presence of the stellar tail to explain the 
anti-correlation between the H i and optical light in the 
N clump. Whatever their cause, these observations point 
to some important differences occurring between gas and 
stars during the formation of some tidal tails. 

4.4. Conditions for Luminous Infrared Phase 

A main objective of our studies on the atomic gas in IR 
luminous galaxies (L/^, > SxIO^^Lq, Hibbard & Yun 1996 
and in preparation) and the molecular gas in the less-IR lu- 
minous systems of the Toomre Sequence (Yun & Hibbard 
1999) is to investigate the conditions necessary for fuel- 
ing a highly elevated period of star formation. Certainly a 
high gas content is one of the pre-requisites for fueling such 
large starbursts (Sanders et al. 1988). We are interested 
in knowing whether all mergers between gas rich galaxies 
pass through such a phase, or if there are requirements 
on the progenitor properties or encounter characteristics 
to trigger such bursts. To gain some insight into how this 
efficient mode of star formation is triggered in Arp 299, we 
draw upon statistical studies of luminous IR galaxies and 
the results of numerical simulations and to interpret our 
observations. 

The tight correlation between molecular gas column 
density {T,h2) and L/^, especially when compared with 
the much looser correlation of gas mass with Lju, shows 
that IR luminous starbursts are always precipitated by a 
corresponding increase in the nuclear column densities of 
molecular gas (Scoville et al. 1994; Solomon et al. 1997; 
Kennicutt 1998; Taniguchi & Ohyama 1998; Yun & Hib- 
bard 1999). The CO observations of Arp 299 (Sargent & 
Scoville 1991, Aalto et al. 1997, Casoli et al. 1999) show 
that the highest column densities in Arp 299 occur at the 
nucleus of IC 694, with Y,h2 over three times higher than at 



the nucleus of NGC 3690. Additionally, HCN mapping ob- 
servations (Aalto et al. 1997, Casoli et al. 1999) shows that 
the vast majority of the very dense gas (n > lO^cm"'^) is 
associated with the nucleus of IC 694. Since the two pro- 
genitor galaxies appear to be of nearly equal mass and of 
similar total gas content (Table P, we suggest that the 
reason there has been so much more gaseous dissipation 
towards the nucleus of IC 694 is due to its retrograde spin 
geometry (see also Sanders et al. 1988). In such encoun- 
ters, the gas within the disk is not pulled into tidal bridges 
and tails, but rather remains in the disk where it feels the 
perturbation from the second system alternately pulling it 
outward and pushing it inward during an orbit (e.g. TT72; 
White 1979; Noguchi 1991). This excites large epicyclic 
motions in the gas, which in turn leads to higher collision 
rates (e.g. Olson & Kwan 1990; Noguchi 1991) and hence 
to more gaseous dissipation. 

This suggestion would seem to be at odds with the re- 
sults from numerical simulations, which find only a mod- 
erate dependence of gaseous dissipation and/or starforma- 
tion activity on the encounter spin geometry (e.g. Mihos & 
Hernquist 1996; Barnes & Hernquist 1996). Further, these 
simulations invariably have the starbursts occurring only 
when the nuclei finally coalesce, while in Arp 299 ( and 
other IR luminous mergers, Murphy et al. 1996, Mihos & 
Bothun 1998), the nuclei are still well separated {d > 4.7 
kpc) . A possible solution to these concerns is found in the 
simulations of Noguchi (1991), who stresses the episodic 
nature of starbursts in mergers (see also Olson & Kwan 
1990). Noguchi notes that while the maximum enhance- 
ment of activity is not strongly dependent on the spin 
geometry of the encounter, the time evolution of the ac- 
tivity is. In particular, his prograde-retrograde encounter 
produces a longer period of enhanced activity over other 
spin combinations and also leads to a slower coalescence of 
the nuclei. As a result, the cloud collision rate is still sig- 
nificantly increased as the nuclei separate (see his Fig. 9). 
This model predicts that the most luminous phase of Arp 
299's starburst is still in the future, in which case this 
merger is predicted to out-shine even Arp 220. 

However, the numerical simulations are necessarily ad- 
hoc in nature, and the micro-physics of hydrodynamics 
(e.g. SPH as in Mihos & Hernquist and Barnes & Hern- 
quist vs. sticky particle as in Noguchi), star formation, and 
"feedback" (energy input back into the interstellar medium 
due to star formation) is poorly understood. Therefore, 
before drawing too many inference from the simulations, it 
would be prudent to first test the numerical formalisms for 
star formation, dissipation, and feedback by direct com- 
parison with real systems. Until simulations can repro- 
duce starbursts occurring at the times and locations as 
observed, we must be cautious when extrapolating their 
results too far into the future. 

The observational fact that systems can experience 
much of their star formation when the nuclei are well sep- 
arated has important ramifications for the structure of the 
resulting remnant. In particular, stars are dissipationless 
and cannot radiate their orbital energy away. Therefore 
stars that form while the nuclei are well separated will be 
spread over a larger range of radii in the resulting rem- 
nant than those that form when the nuclei are practically 
merged. This population can experience violent relaxation 
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(e.g. Barnes 1992; Barnes & Hernquist 1996), resulting in 
an extended dynamically hot population. The resulting 
remnant will have a much more regular luminosity profile 
than a merger in which the gas all sinks to the very center 
of the potential before forming stars (c/. concerns raised 
in Mihos & Hernquist 1994). 

4.5. Lack of Tidal Dwarfs 

The widespread lack of small-scale correlation between 
optical and gaseous density enhancements within the tidal 
tail was unexpected. In previous H i and optical observa- 
tions of six double-tailed merging systems (van der Hulst 
1979; Hibbard & van Gorkom 1993, 1996; Hibbard et al. 
1994), five dwarf galaxy like concentrations of gas, light 
and Ha emission appear entrained within tidal tails and 
coincident with distinct features in the optical color maps 
and the H i velocity dispersion maps. These concentra- 
tions have been interpreted as "tidal dwarf galaxies" form- 
ing out of the expanding tidal material (see also Schweitzer 
1978; Mirabel, Dottori & Lutz 1992; Due et al. 1997; 
Malphrus et al. 1997), an idea that has been supported 
by numerical simulations (Barnes & Hernquist 1992, 1996; 
Elmegreen, Kaufmann & Thomasson 1992; Hibbard & Mi- 
hos 1995). In addition to these dwarf-sized concentrations, 
the previously observed tidal tails often contain numer- 
ous knots of gas, light, and star-forming regions (see also 
Hutchings 1996; Hunsberger, Charlton & Zaritski 1996; 
Mihos & Bothun 1998). We therefore fully expected sim- 
ilar features to appear in a tidal feature as long and as 
massive as that seen in Arp 299. 

However, the optical tail of Arp 299 is remarkably 
smooth and featureless (Fig. ||). There are a number of 
point sources falling upon the optical tail, but the surface 
density of such points is similar to that of point source 
elsewhere in the field. Their colors are generally quite red 
{B — R > 1.7, Fig.j^a), and most are probably background 
galaxies. And while there is considerable structure within 
the H I tail, there are no H i peaks with a corresponding en- 
hancement in the underlying starlight (Fig. ^). Further, 
there are only three regions which exhibits a significant 
(>20%) increase in the Hi velocity linewidth (indicated 
in Fig. 0d) — along the southern edge of the N clu mp and 
at the Hi knot in the inner tidal filament (||3.2) — and 
these features do not correspond to any notable features 
in the optical images or color maps. 

To further investigate this question, we estimate the im- 
portance of self-gravity within the tail by comparing the 
luminous mass (M/„m) to the dynamical mass (Mdyn)- 
The luminous mass is estimated by adding the gas mass 
and the stellar mass. The latter is estimated from the 
B-band surface brightness and adopting a stellar mass to 
light ratio of {M/L)b = 2MqLq^, characteristic of stel- 
lar disks (Bottema 1993). The dynamical mass is calcu- 
lated via Mdyn — 1-76 x 10^ Af© x a'^j x ri/2 (Binney 
& Tremain 1987; see also Kaufman et al. 1997), where 
ri/2 is the FWHM radius of the beam (~ 1.9 kpc for the 
high resolution data). This will over-estimate the effects 
of self-gravity if (1) the line-of-sight dimension of the tail 
is larger than the beam width (~ 3.8 kpc), in which case 
the above prescription overestimates the stellar mass; and 
if (2) the stars contribute significantly to Mium and have 

^^Using a velocity dispersion of 20 km s^^ as inferred for the stars 



a velocity dispersion greater than aHi (e.g. van der Kruit 
1988). Since we suspect both of these conditions to hold, 
our calculation of Mium/Mdyn should be a stringent upper- 
limit. We derive a lower limit to Mium/Mdyn by setting 

Using the high resolution data, we find that 
Mium/Mdyn < 0.3 everywhere, with the highest val- 
ues on the brighter optical peaks of the inner tidal fil- 
ament^^. The peaks in the Hi column density along 
the outer filament and in the N clump have 0.2 < 
Mium/Mdyn < 0.3, and the Hi knot in the inner fila- 
ment has 0.1 < Mium/Mdyn < 0.2. For comparison, 
the tidal dwarfs in the tail of NGC 7252 (Hibbard et al. 
1994), NGC 3921 (Hibbard & van Gorkom 1993, HvG96) 
and NGC 4038/9 (Mirabel, Dottori & Lutz 1992) have 
0.50 <Mium/Mdyn ^1- Therefore Arp 299 lacks the po- 
tentially self-gravitating dwarf-like concentrations seen in 
other tailed systems, and we conclude that tidal dwarf for- 
mation is not a ubiquitous process. 

This suggests that the formation of such sub-structures 
may depend on factors other than the ability to raise a 
decent sized tail. We note that other tidal tails which 
contain candidate tidal dwarf galaxies have similar gas 
column densities as in the tail of Arp 299, but have op- 
tical tails that are ~2 mag arcsec"^ brighter. Addition- 
ally, in these systems the gas and optical material coin- 
cide. Both of these factors increasing the underlying mass 
density within the tail. Further, since the gas in observed 
tidal dwarfs are a significant contribution to the dynamical 
mass {Mgas/Mdyn ~ 0.5), it may be the case that what- 
ever proce ss removes the H i from the inner tidal filament 
(see § 4.3.2| ) also suspends the dwarf formation process. 

Whatever the explanation, the present observations sug- 
gest that one should be cautions about simply measuring 
the mass or light contained within the highest gas column 
density contours or optical isophotes and calling these con- 
centrations distinct entities. 



5. CONCLUSIONS 

• A new tidal tail has been discovered in the on-going 
merger Arp 299. The tail stretches 180 kpc to the 
north, and must have taken at least 750 Myr to 
form. This is much longer than the age of the current 
star burst, suggesting that the starburst and merger 
clocks start at different times and run at different 
rates. If the present starburst began when the tidal 
tail was first launched, the spread of ages in the re- 
sulting remnant will be much larger than is usually 
assumed. 

• What had previously been interpreted as a single 
rotating gaseous disk surrounding the two merging 
systems is shown to be a gaseous disk associated with 
the eastern system, IC 694. The kinematic informa- 
tion in concert with the optical morphology allows 
us to deduce that NGC 3690 experienced a prograde 
encounter while IC 694 experienced a retrograde en- 
counter. From the cold gas contents and NIR mor- 
phologies, we suggest late type progenitors for both 
systems (Sab-Sb for IC 694, Sbc-Sc for NGC 3690). 



in S3.2 would give Mi^^/M^y„ < 0.04 on this feature 
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• The observations show that Arp 299 is experiencing 
an IR bright starbursts while the nuclei are still well 
separated (dr^4.7 kpc). We suggest that this is due 
in part to the retrograde spin geometry of IC 694. 

• The tidal morphology is unique among mergers: 
there is a low surface brightness stellar tail with 
low column density H i and a conspicuous parallel 
gaseous tail of higher H i column density and with 
little, if any, associated starlight. The tails have very 
similar kinematics suggesting that they are physi- 
cally connected. We suggest that such bifurcated or 
parallel tails result from a progenitor with a warped 
disk. However, this explanation does not explain the 
low H I content of the optical tail. 

• The similarity between this morphology and that of 
other starburst galaxies with direct evidence of su- 
perwind outflows (M82, NGC 4631, NGC 520, Arp 
220) suggest that the galactic superwind recently im- 
aged by Heckman et al. (1998) may be responsible for 
clearing the optical tail of much of its cold diffuse gas. 
Local ionizing sources within the tidal tails may also 
play a role in some of the observed optical/gaseous 
anti-correlation. These issues are explored in more 
detail in HVY99 

• The absence of dwarf galaxy-like mass concentra- 
tions in the tidal tail suggests that tidal dwarf forma- 
tion is not a ubiquitous process. It is possible that 
the whatever process removed gas from the optical 
tail may have consequently suppressed the dwarf for- 
mation process. 
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APPENDIX 
WHO IS IC 694? 

There is some confusion in the printed and electronic literature as to which systems actually comprise Arp 299 (e.g. lAU 
circular 6859; NED)^. Arp does not offer many clues, and even confuses the matter by incorrectly designating Arp 296 
as NGC 3690+IC 694 in Tables 1 & 2 of his Atlas, while giving Arp 299 no further designation (Arp 1966). However the 
description in Arp's Table 1, and the order of his flgures (Arp 296 appearing next to other bridge-tail systems) makes it 
clear that Arp published the pictures under the correct number, and gave them the correct description in Table 1, and 
just got the "Designation" switched in both tables. 

The description given for NGC 3690 in the NGC (Dreyer 1888) translates as "pretty bright, pretty small, very little 
extended at PA=80, pretty gradually brighter towards the middle, small stars south following near" . According to the 
NGC "pretty small" means a diameter of 50"-60", and this radius would incorporate much of the brighter parts of the 
system (including the nucleus of the eastern system). Given that the five highest surface brightness objects are within 
7" of source B in the western system, through a small telescope this object must have looked like a single nucleus at 
the location of source B in Fig. ^ with a cloud of surrounding nebulosity, and it is unlikely that two nuclei could have 
been discerned. It therefore seems likely that the original designation refers to the entire Arp 299 system, although only 
one nucleus was identified (corresponding to the western system, what we have been calling NGC 3690 in the above). 
Therefore whether NGC 3690 should strictly apply to the entire system or to only the western system depends on if you 
believe galaxies are identified by their nuclei or their integrated light. We remain agnostic on this point, and from here 
investigate which object was designated as IC 694. 

The earliest and clearest association of IC 694 with the eastern system can be traced to Vorontsov-Veryaminov, in his 
Atlas and Catalogue of Interacting Galaxies (Vorontsov-Vel'yaminov 1959), where Arp 299 appears as VV118a-e. Source 
VV118a is the galaxy to the East (what we have been referring to as IC 694), VVllSb the system to the west (what we 
have been calling NGC 3690), VVllSc is MCG+10-17-2a, the compact spheroidal 1' to the NW, and VV118d & e are 
bright HII regions in the outer disk regions to the NW. In the accompanying table to the Atlas, VVllSa is identified 
as IC 694, VVllSb is identified as NGC 3690, and sources c,d or e have no other designations. The vast majority of 

^''For this discussion, it is instructive to look at an earlier photographic representation of the Arp 299. For this, we refer to reader to Plate 1 
of Morgan (1958). In this figure, the NW spheroidal appears as a very faint fuzz well separated from and much dimmer than either of the two 
interacting disks, which appear as distinct objects. 
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researchers since have foUowed this designation. Nilson first brings up the question of the identification of fC 694 with 
VV 118a in the UGC (Nilson 1973), in his notes on UGC 6471/2: "Arp 296 (identified as NGC 3690 + IC 694 by Arp) 
are faint objects north-preceding UGC 06471 + UGC 06472, magnitude approximately 17 and 21 identification of IC 694 
uncertain, may be a small object north-preceding and inside the outer parts of the double system" . [The use of the name 
Arp 296 is immediately attributed to Arp's error in making his tables, discussed above]. 

IC 694 was discovered by Swift (1893), whose discovery description translates as "close double with NGC 3690 — object 
247 of list I of Sir William Herschel. Suspected at 132x magnification, verified at 200x". In the Index Catalogue of 
Nebulae (IC, Dreyer 1895) this becomes "very small, forms double nebula with no. 247 of list I of Sir William Herschel". 
"very small" translates as 10"-20" diameter, which is on the large side for the NW spheroidal {0fwhm^4:")- There is no 
definition of the term "double" in the IC^, and perhaps this is where the confusion lies. Certainly with todays knowledge 
we would call Arp 299E-I-W a double, and if anything the spheroidal to the NW as a possible companion. However, is 
this how Dreyer and Swift used the term? 

So the question is, when Swift identified IC 694, did he subdivide the NGC 3690 entry, or add an entry for the 
faint nebulosity to the NW? The position for IC 694 is given as 11:20:44 -^59:20 (epoch 1860) in the IC, while the NGC 
position for NGC 3690 is 11:20:45 +59:19, i.e. IC 694 lies approximately 1' to the NW of NGC 3690. This seems like pretty 
conclusive evidence that IC 694 indeed refers to the spheroidal to the NW. However, Dreyer remarks in the introduction to 
the IC that Swifts positions are "generally reliable within one or two minutes of arc, but larger errors occur occasionally..." 
We also learn from Dreyer that Swifts observations were taken through a 16" refractor at the Rochester Observatory. 
Since the NW galaxy is 3 mag fainter than either of the disk galaxies (tob=16.2 mag from CCD observations reported 
above), is it even possible for Swift to have seen the NW companion? If not, it is likely that he simply resolved NGC 
3690 into two separate condensations of nearly equal luminosity, adding a separate catalogue entry for the second system. 
A similar super-classification appears to have happened in the case of NGC 4861/IC 3961 (Arp 266) and NGC 2207/IC 
2163. 

Recent attempts to repeat these observations with similar sized telescopes by amateur astronomers suggest that (1) 
the second light concentration now associate with the nucleus of Arp 299E is very easy to discern as distinct from Arp 
299W; and (2) the NW spheroidal is indeed a very difficult, but not an impossible, object. Combined with the correct 
relative position of IC 694 given by Swift, it seems that there is strong evidence that IC 694 properly designates the 
NW compact spheroidal. However, since these observations were made with prior knowledge of both the existence and 
location of the NW spheroidal, this confirmation is not entirely independent. Since it is nearly impossible to prove this 
issue conclusively one way or the other, we choose to follow the vast majority of researchers and use the designation 
introduced by Vorontsov-Vel'yaniinov (1959), concluding as Sulentic & Tift (1973) that "Any further question as to which 
objects Dreyer was referring to can only be of historical interest" . 

JEH offers special thanks to Harold Corwin, Steve Gottleib, Malcom Thomson, Bill Vacca and Dennis Webb for many 
useful discussions on this issue. 
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Fig. 1. — B-band image of the Arp 299 system. The optical image is displayed with a logarithmic transfer function, with 
white set to 28.5 mag arcsec^^ and black set to 22.75 mag arcsec"^. IC 694 is to the left (east) and NGC 3690 is to the 
right, and north is up. Regions denoted as A,B, and C-C by Gehrz, Sramek & Weedman (1983) are so labelled. We also 
identify the associated system MCG+10-17-2a and the unrelated background bar/bridge/tail system Arp 296. A scale bar 
is drawn indicating 10 kpc, assuming a distance of 48 Mpc to Arp 299. The putative nuclei of the progenitors, as judge by 
NIR imaging and spectroscopy, are at regions A and B. Both of these locations are the sites of very dense molecular gas 
concentrations. The label C-C indicates the region of disk overlap, which is also the site of a significant concentration of 
molecular gas. 
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Fig. 2. — Position- Velocity profiles of the inner regions of Arp 299, demonstrating the kinematics of the inner H i features. 
The top left panel presents the i?-band image displayed as in Fig. ^, with contours from the high resolution H i integrated 
intensity map superimposed (FWHM=17"xl5", as indicated by the ellipse in the lower left corner of the panel). Negative 
contours are drawn as dotted white lines, indicating the region of H i absorption. The positive contours begin at 11.05 mJy 
beam~^ km (5x lO^'^^cm"^), while negative contours begin at -33.15 mJy beam~^ km s~^. Successive contour levels 
are twice the previous level. The black dotted contour indicates a column density of IxlO^^cm"^ from the low-resolution 
Hi data {Ofwhm—S^"). To the right is a position- velocity map of the high resolution data after summing the emission 
in right ascension, and below is a similar plot, made by summing the emission in declination. Positions are measured 
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Fig. 3. — H I channel maps of the inner regions of Arp 299 contoured upon the _B-band image. The higher resolution 
data are used after Hanning smoothing in velocity by a factor of two to a channel spacing of 21.4 km s~^. The beam 
size (17"xl5") is indicated in the lower left hand corner of the first panel, and each panel is labelled with the channel 
number and heliocentric velocity. Dotted white lines indicate negative contours, which are drawn at lax [-3, -5, -7, -10], 
where 1(7=0.25 mJy beam~^ is the single channel noise level. Positive contours are drawn at lax [3, 5, 7, 10], with 
higher contours drawn a factor of 1.5 times the previous level. The lowest contour corresponds to a column density of 
7.2x101^-2^ 
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Fig. 4. — Position- Velocity profiles of the entire Arp 299 system, demonstrating the kinematics of the tidal H i features. 
The top left panel presents a false color and contoured image of the intermediate resolution data cube (22"x20"). Features 
referred to in the text are labelled. The adjacent panels are position-velocity profiles, constructed by summing the data 
cube in right ascension (right) or declination (bottom). Values below la were not not used when constructing the position- 
velocity plots. Contours are drawn starting at a level of 8 mJy bcam^^ km s^^ (2x lO^^cm^^), with higher contours drawn 
a factor of 2 times the previous level. 
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Fig. 5. — Deep optical image of the entire Arp 299 CCD mosaic, showing the total extent of the stellar tail. The deep B- 
and _R-band data arc shown on the left and right, respectively, using a logarithmic transfer function. The upper panels 
(a &: b) show a greyscaled representation of the data, while the lower panels add contours. The greyscales give a better 
feeling for the structure at the lowest light levels, while the contoured images provide a quantitative measure of the light 
peaks. The greyscales cover the range /Ub=28.5 mag arcsec"^- 22.75 mag arcsec"^ and ij,r=28.5 mag arcsec"^- 22.0 
mag arcsec"^. In the upper panels, white contours are drawn every 1 mag arcsec"^ starting at /iB=24 mag arcsec"^ and 
fiji—25 mag arcscc^^. In the lower panels black contours are added every 1 mag arcsec"^ starting at /is=28 mag arcsec"^ 
and hr=27 mag arcsec"^. Three stars referred to in the text are labelled in panel (a). 



Fig. 6. — An RGB representation of the Arp 299 optical and HI data. Yellow and green indicates optical emission, while 
blue indicates H I emission. Regions containing both forms of emission appear as the sum of these colors. This figure 
emphasizes the differences between the stellar and gaseous tidal morphologies. 
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Fig. 7. — Detail of the Arp 299 tidal tail, showing the relationship between the gaseous and stellar morphologies and 
kinematics, (a) B-R color map constructed from the smoothed data using only pixels with a signal to noise greater than 
or equal to 3 in both the B- and _R-band images. A spectral transfer function is used with the following color mappings: 
Purple: B-R <-0.2 mag, blues: B-R = -0.2 - +1.0; green-yellow: B - i? = 1.0 - 2.0 mag; red: B-R >2.0 mag. (b) 
False color rendition of the _B-band data with contours of the intermediate resolution data overlaid. Contours are drawn 
at levels of 4 mJy km s~^ beam"'^ (IxlQ^^cm"^) x[3, 6 10, 15, 20, 25]. (c) intensity weighted mean Hi line-of-sight 
velocity. The intermediate resolution has again been used, but with a higher threshold for the moment algorithm. The 
color boundaries are indicated with labels giving the velocity (in km s^^) with respect to NGC 3690 {Vhei «3040 km s~^, 
§ ^■41 ). Contours are the same as in (b). (d) Hi line-of-sight velocity dispersion, using the same data and threshold as in 
(c). The colors correspond to the following velocity dispersion ranges: purple-blue = 5-8 km s~^; green = 8-10 km s~^; 
yellow-orange = 10-13 km s~^; red = 13-20 km s~^. The regions of high velocity dispersion discussed in §3.3 and §4.^ 
are indicated. 
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Fig. 8. — Channel maps of the tidal regions of Arp 299 contoured upon the S-band image. The grayscales are drawn 
from /iB=28.5 mag arcsec^^ (white) to 25 mag arcsec"^ (black). The more sensitive low-resolution Hi cube is used at 
full velocity resolution (10.52 km s~^ channel width). The 35" beam size is indicated in the lower left hand corner of 
the first panel, and each panel is labelled with the channel number and heliocentric velocity. Dotted yellow lines indicate 
negative contours, which are drawn at Ictx [-3, -5, -7, -10], where 1(7=0.37 mJy beam~^ is the single channel noise level. 
Positive contours are drawn at levels of lax [3, 5, 7, 10, 15, 22.5]. The lowest contour corresponds to a column density 
of IxlO^^cm^^. Only channels showing Hi emission from the tail are shown. 
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Fig. 9. — Position- Velocity slices taken along both gaseous filaments. The left panels show the location of the slice on the 
low-resolution H i iintegrated intensity maps, while the right panels show the corresponding position- velocity profile. For 
the latter, distance along the slice is plotted along the y— axis, and velocity with respect to 3100 km is plotted along 
the x— axis. Crosses are drawn at intervals of 100" along the slice in the left panel, and arrows indicate the corresponding 
locations in the right panel. Specific regions referred to in the text are labelled, (a) & (b) Position- Velocity slice taken 
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Fig. 10. — Greyscale image of the companion to Arp 299, taken from Version II of the Digitized Sky Survey, with Hi 
contours superimposed. The higher resolution data is used with the 17"xl5" beam indicated by the ensquared ellipse to 
the lower left. Hi column density contours are drawn at [5, 10, 15, 20, 30] xlO^^cm"^. 



Fig. 11. — H I channel maps of the companion to Arp 299 contoured upon a greyscale of Version II of the Digitized Sky 
Survey. The higher resolution (17"xl5"), full velocity resolution (Au=10.52 km s""'^) Hi data are used, with the spatial 
resolution indicated by the ensquared ellipse to the lower left of the first panel. Contours are drawn at la x[-3, 3, 5, 7], 
where 1(t=0.20 mJy beam~^ is the single channel noise level, corresponding to a column density of 9.2x lO^^cm"^. 
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Table 1 
VLA Observing Parameters 



Parameter 



ll'»25'" 44!2 +58°50'18" 
3080 km 

30' 

1203+645 
3C268 
3.125 MHz 
63 

10.5 km s"^ 



Array Configuration 


D 


C 


C+D 


Date 


4/18/95 


12/10/94 




Synthesized Beam 








— FWHM: Major Axis x Minor Axis 


56"x46" 


17"xl5" 


22"x20" 


— Position Angle (East of North) 


+37° 


-22° 


-12° 


Time on Source (hrs) 


2.1 


2.9 


5.0 


Noise Level (la) 








— Flux Density (mJy beam~^) 


0.56 


0.37 


0.30 


— Column Density (xlO^^cm"^ beam~^ ch~^) 


0.25 


1.7 


0.80 


— Brightness Temperature (K beam~^) 


0.13 


0.88 


0.41 



Table 2 

Optical Observing Parameters 



Parameter 



Telescope 




UH 88" 


Detector 




Tek 2048 


Readout Mode 




binned 1x1 


Focal Ratio 




f/10 


Pixel Size 




0'.'22 


Field of View: 






— Single CCD frame 




7.'5 X 7'.5 


— Final Image 




8.'7 X 12'A 


Date 


6/03/95 


01/07/96 


Filters 


R 


B 


Effective Exposure Time 


~ 2 X 600s 


~ 2 X 900 s 


Seeing 


1'.'2 


o'.'g 


Sky Brightness (mag arcsec"-^) 


20.9 


22.9 


Ict Sky Noise^ (mag arcsec"^) 


25.8 


27.0 


30- Sky Noise'' (mag arcsec"^) 


26.9 


28.0 



^la noise in original unbinned image 

^3(7 noise after binning 9x9; this should be close to the 3a detection limit, and we can often trace extended features to 
0.5 mag fainter than this. 



Phase Center (aigso (^1950) 
Velocity Center (Heliocentric) 
Primary Beam (FWHM) 
Phase Calibrator 
Flux Calibrator 
Bandwidth 
Number of Channels 
Channel Separation 
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Table 3 

Measured Properties of the Arp 299 System 



Quantity 



Units 



Total 



Disk* 



Tail 



Companion'' 



Optical: 

Lb' 
Lr'= 
B 
R 

H I, measured: 
Velocity Range** 
J Sdv 
Mh^ 
Mhi/Lb 
Mhi/Lr 

H I, corrected? 

/ Sdv 
Mhi 
Mhi/Lb 
Mhi/Lr 



{Lq,b) 
{Lq,r) 
(mag) 
(mag) 



(km s-i) 
(Jy km s"-') 
{Me) 
{MeL^') 
{MeL^') 



(Jy km s ^) 

{MqL-') 
{MqL-') 



4.5 X IQi" 
4.9 X IQi" 
12.31 
11.04 



2800 



3380 
17.7 
9.6 X 10^ 
0.2 
0.2 



21.7 
1.18 X 10^° 
0.3 
0.3 



4.3 X 10^° 
4.7 X 10^° 



2910 



3350 
11.7 
6.3 X 10^ 
0.15 
0.13 



15.7 
8.5 X 10^ 
0.20 
0.18 



1.6 X 10^ 
1.2 X lO'^ 



3020 



3210 
6.1 

3.3 X 10^ 
1.6 
1.7 



2 X 10^ 
2 X 10^ 
18 
17 



3230 



3250 
0.14 
8 X 10^ 
0.04 
0.04 



"^Disk values include contribution from minor axis clump and western plume 

^HI detected companion located at aigso = ll''24™28?0, i5i95o = +58°48'05" and Whei = 3245 km s'^ 

■^Luminosities calculated assuming ^=+5.48 and M0,fl=+4.31. Values for disk are measured to /iB=26mag 

arcscc^^, /i/j=25.5mag arcscc^^. Values for tail arc measured after removing point sources. The values for the com- 
panion are estimated from the DSS image and arc very approximate; see text. 

■* Range of H i velocities (heliocentric). Taken from the first- moment image except for the total, which is the range seen 
in absorption. The uncertainty is ±5 km s~^. 

"Integrated Hi mass, calculated using Mhi = 2.356 x lO'''Af0 J Sdv, where A is the distance in Mpc, and J Sdv is 
the integrated H i emissivity, in Jy km s~^. Following Sanders, Scoville & Soifer (1991), we adopt 48 Mpc as the distance 
to Arp 299. 

^Absorption corrected value; see text. 
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Table 4 
Global Quantities 



Quantity 


Units 


Total 


IC 694 


NGC 3690 


Lb 


Lq,b 


4.5 X 10^° 


2.8 X 10^° 


1.7 X 10^° 


Mhi 


Mq 


1.18 X 10^° 


6 X 10^ 


6 X 10^ 


Mh,^ 


Mq 


1.7 X 10^° 


1.1 X 10^° 


6 X 10^ 


Mgas"" 


Mq 


3.3 X 10^° 


1.9 X 10^° 


1.4 X 101° 


Mhi /Lb 


MqL-^ 

Mq Lq^b 


0.3 


0.2 


0.4 


MhJLb 


0.4 


0.4 


0.4 


Mgas/LB 


MqL-]j, 


0.7 


0.7 


0.8 


MhJMhi 


1.4 


1.8 


1.0 


Lir'^ 


(Lqmi) 


7.9 X 10" 






Lir/ Lb 




18 






Lir/Mh, 


{LqM^') 


69 







^From Casoli ct al. (1999), converting to our distance and a CO-to-i/2 conversion factor of 3 x 10^° cm^^ (Kknis^i)^^. 
We have equally divided the gas located at the overlap region C — C between IC 694 and NGC 3690, and assigned the 
"widespread" emission to IC 694. 

^Total gas mass corrected for the presence of He: Mgas = Mh2 + 1-34 x Mhi 



^Far infrared luminosity over the range 8 — lOOO/xm, from Sanders et al. (1991). 
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